Investigation of X-ray cavities in the cooling flow system Abell 1991 
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ABSTRACT 



We present results based on the systematic analysis of Chandra archive data on 
the X-ray bright Abell Richness class-I type cluster Abell 1991 with an objective to 
investigate properties of the X-ray cavities hosted by this system. The unsharp masked 
image as well as 2-d (3 model subtracted residual image of Abell 1991 reveals a pair of 
X-ray cavities and a region of excess emission in the central ~12 kpc region. Both the 
cavities are of ellipsoidal shape and exhibit an order of magnitude deficiency in the X- 
ray surface brightness compared to that in the undisturbed regions. Spectral analysis 
of X-ray photons extracted from the cavities lead to the temperature values equal to 
1.77i£^ keV for N-cavity and 1.53±g;g| keV for S-cavity, while that for the excess X- 
ray emission region is found to be equal to 2.061qq7 keV. Radial temperature profile 
derived for Abell 1991 reveals a positive temperature gradient, reaching to a maximum 
of 2.63 keV at ~ 76 kpc and then declines in outward direction. 0.5—2.0 keV soft band 
image of the central 15"region of Abell 1991 reveals relatively cooler three different 
knot like features that are about 10"off the X-ray peak of the cluster. Total power of 
the cavities is found to be equal to ~ 8.64 x 10 43 erg s -1 , while the X-ray luminosity 
within the cooling radius is found to be 6.04 xlO 43 erg s _1 , comparison of which 
imply that the mechanical energy released by the central AGN outburst is sufficient to 
balance the radative loss. 

Subject headings: galaxies:active; galaxiesxlusters; X-rays:galaxies:clusters; cooling flows: 
intergalactic medium 



-3- 



1. Introduction 



High-resolution X-ray observations of bright cool-core galaxy clusters with Chandra and 
XMM-Newton telescopes have shown that, even in the cores of clusters with high cooling ra te, 



the observed quantum of gas cooled below T < 1 keV is significantly small (IZuHone 



20111) . 



Further, the rate at which mass is being deposited in the cores of the dominant cooling flow 
clusters is found to be an order of magnitude sm aller than that estimated using standard models 



(Peterson et al 



2003 



Peterson and Fabian 



2006|) . This mea ns some compensatory steady heating 



mechanism must be operative in the cores of such clusters (|McNamara et al. 



2006). A variety of 



mechanisms have b een proposed to compensa te the radiative losses, that include, the magnetic 



field re-connection (Soker anc 



(Naravan and Medvedev 



Sarazin 



(| 19901) ). thermal conduction due to electron collisions 



200 ll). turbulent conductio n (|Voigt and Fabian 



cosmic rays (|Colafrancesco and Marchegiani 



2004|) and heating by 



20081) . However, the most promising source of 



heating for replenishing the radiative loss of the hot gas in clusters and groups of galaxies is the 
energy input from the super-massive black hole (SMBH) that is residing at the core of such c luster 



(IBoehringer et al. 



19931: 



Binnev and Tabor 



19951 DVIcNamara et al. 



2006; 



Fabian et al. 



2006). 



One of the early clues explaining role of the active galactic nucleus (AGN) in reheating of the 
inter-galactic medium (IGM) was the high detection rate (~70%) of radio activity in the central 



galaxy of the X-ray bright clusters relative to the X-ray faint clusters (|Burns 



1990) 



X-ray deficient cavities detected in numerous groups and clusters of galaxies provided with 



the direct evidence for the AGr> 



Dong et al. 



20101 : 



Randall et al. 



feedb ack in such systems (jBirzan et al. 



2004; 



David et al. 



2009; 



201 lb . Often such cavities are seen near the center of clusters in 

buoyantly rising bu bbles in 



about few to few tens of kpc aro und the X-ray center and appear as 



the intra-cluster medium (ICM) (|Dong et al. 



20101 : 



David et al. 



2009 



Pandge et al. 



such cavities are found in association wi th radio jets or bubbles / lobe s and are often fi 



20121) . Further, 
led with 



relativistic particles and magnetic fields (|McNamara and Nulsen 



2007 



Gitti et al. 



2012 ). Wealth 
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of the observational evidences over last few decades have shown that, roughly 70 - 75% of the 
X-ray bright cool-core clusters harbor detectable number of cavities with energy input from the 



centra l AGN roughly sufficient to compensate the cooling of gas in the clusters (|Rafferty et al 



20061) . This means, feedback from the central AGN plays an important role in shaping the 
morphology of the hot gas and also in the thermodynamical evolution of the galaxy clusters 



( Guo and Oh 



20081) . providi ng with a reliable means to estimate the mechanical energy injected 



by the SMBH in to the ICM (IJones and De Young 



2005). 



In this paper we present results based on the analysis of Chandra observations of the 
X-ray bright cluster Abell 1991 with an objective to investigate properties o f the X-ray cavities 
hoste d by this system. Abell 1991 is classified as the Bautz-Morgan type I (|Bautz and Morgan 
19701) c luster with a centrally do minant galaxy and is defined as the Abell Richness Class I type 



system ( 



Struble and Rood 



1999J). Structure of the paper is as follows: Section 2 describes X-ray 
observations and data preparation methodology. Section 3 discusses the imaging analysis with 
an emphasis on the investigation of X-ray cavities. In Section 4 we discuss spectral properties of 
hot gas within X-ray cavities and other regions of interest, while Section 5 discusses energetics 
associated with the central engine and correspondence between the X-ray and radio observations. 
Summary of the results from the study is presented in Section 6. We assume H =73 km s _1 Mpc -1 
in this paper. At the red-shift of Abell 1991, luminosity distance is 255 Mpc and corresponding 
angular size is 1.17 kpc arcsec" 1 . 



2. Observations and Data Preparation 



Abell 1991 was observed with Chandr a Advanced CCD Im aging Spectrometer (ACIS-S) 



on 16-17 December 2002 (ObsID 3193, see 



Sharma et al 



(|2004r) for details) for a total exposure 



of 38.8 ks. Event-2 data sets available in the archive of Chandra observatory were acquired 
for the present study and were reduced following the standard tasks available within the X-ray 



analysis package CIAO V 3.4.cQ in conjunction with the Chandra calibration data base (CALDB) 
V 3.4.0 provided by the Chandra X-ray CenteJ^. Initially the event files were screened for 
cosmic rays using ASCA grades and were reprocessed to apply most up-to-date corrections for 
the time-dependent gain change, charge transfer efficiency, and degraded quantum efficiency of 



the ACIS detector. Though partic 



have been detected (ISharma et al. 



e bac kground during observations was stable and no flares 



2004|) . however, as a routine process we performed light 



curve filtering to detect the time intervals affected by the background counts exceeding 20% 

of their mean count rate. This resulted in the net effective exposure time of 38.3 ks. X-ray 

J] 

background component was adequately modeled using the "blank-sky" data setqj, that included 
both particle-induced and unresolved sky components. Spectra from regions of interest during 
this study were extracted using the AC IS SPEC task and corresponding Redistribution Matrix 
Files (RMF), Ancillary Response Files (ARF) were generated. The point source detection was 
performed on CCD ID 7 using CIAO WAVDETECT task with the source detection threshold set to 
10~ 6 . This enabled us to detect a total of 21 point sources within S3 chip . We generated exposure 
corrected, background subtracted 0.5 - 3.0 keV energy band image using DMMERGE task available 
within CIAO 3.4.0 and is shown in Fig.[2|a). This figure reveals the bright X-ray region with its 
peak emission centroid located at (RA = 14:54:31.7) and (DEC = +18:38:42), showing an offset 
of ~ 10"rel ative to that in optical/ra dio image (RA = 14:54:31.5; DEC= +18:38:32) centered on 



NGC 5778 (Struble and Rood 



1999) 



'Chandra Interactive Analysis Observation package 

2 http : //cxc . harvard . edu/ 

3 http : //cxc . harvard . edu/contrib/maxim/aci sbg 
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3. X-ray Imaging Analysis 



3.1. X-ray Cavity Detection 



A pair of X-ray cavities in Abell 1991 has already been reported by [Dong et al.1 (|2010|) . 
However, with an objective to enhance visualization and to investigate properties of the X-ray 
cavities in this system, we employed a variety of image processing techniques. These includes, 
deriving unsharp masked image and residual image of Abell 1991. Figure |2](b) shows 0.5-3.0 keV 
Chandra soft band unsharp masked image of Abell 1991 derived after subtracting a 5a wider 
Gaussian k e rnel smoothed image from that smoothed with a narrow 3 a Gaussian kernel (see 

(120121) for details). Resulting unsharp masked image reveals 



Dong et al. 



(2010); 



Pandge et al. 



a pair of X-ray cavities, one on the Northern side (N-cavity) and other on the Southern side 
(S-cavity) of the X-ray center of Abell 1991. Both these cavities are roughly of ellipsoidal shape 
and are located at a projected distance of ~ 12.0 kpc from the X-ray center of the Abell 1991. The 
BCG appears to coincide with the S-cavity and shows an off-set of 10"relative to the X-ray peak. 
As S-cavity has associated radio source, therefore, is defined as "cle ar cavity", while N-cav ity 



failing to show such an association, can be defined as "ghost cavity" (|McNamara et al. 



200 lj). In 



addition to the pair of prominent cavities, an excess X-ray emission is also evident on the Southern 
side of the S-cavity. 

Features evident in unsharp masked image were further confirmed using residual map 
produced after subtracting 2-d smooth model of the Abell 1991 from its original image. A 2-d 
smooth model of Abell 1991 was generated by fitting (3 model to the clean, background subtracted 
image using fitting function Sherpa available within CIAO. Fitting parameters i.e., ellipticity, 
position angle, normalization angle, and local background, etc. were kept free during this process. 
The best fit 2-d model was then subtracted from the cleaned image of Abell 1991 to produce 
its residual map and is shown in Figure |2](c). This figure confirms a pair of ellipsoidal cavities 
evident in the unsharp masked image and exhibit surface brightness fluctuations more clearly. 
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3.2. Azimuthally Averaged Surface Brightness Profile 

Unsharp masked image as well as the 2-d (3 model subtracted residual image revealed a pair 
of X-ray cavities and an excess X-ray emission region (Fig. |2]b & c). To examine the extent of 
flux variation in these regions, we derived surface brightness profile of Abell 1991 using Chandra 
ACIS-S3 data. For this an azimuthally averaged surface brightness profile was computed over 
0.5-5.0 keV range by extracting X-ray photons form concentric circular annuli each of width equal 
to 2.5"starting from 2"up to 175". To produce profile of the "undisturbed" ICM, we excluded the 
cavity regions during this process. The resultant surface brightness profile was fitted with a single 



1986), 



/3-model using the tool SHERPA with the x 2 statistics of Gehrels variance (|Gehrels 

2 -| -3/3+0.5 

: S 1 + I — 

V r o 

where S(r) is the X-ray brightness at the projected distance r, S is the central surface 
brightness, r is the core radius of the X-ray emission and (3 represent the slope of the surface 
brightness profile. Like in other cooling flow systems, distribution of the ICM in Abell 1991 
shows an excess emission in the central region relative to the best fitted model (Fig.|3]/e/? panel), 
whose parameters are r c = 14.50 ± 0.15 kpc and (3 = 0.44 ± 0.01. 

To examine the extent of fluctuations in the surface brightness distribution due to the presence 
of cavities, we considered two conical sectors along the cavity directions i.e., between 0-140° and 
220-360° for N and S cavities, respectively, and extracted X-ray photons from 40 different annuli 
along these sectors. Plots of the count rates along these sectors as a function of radial distance are 
presented in Figure ®(left panel).The profiles derived for both the sectors reveal a significant drop 
in the surface brightness distribution between 8"and 15", roughly of the cavity sizes, compared 
with the best-fit model. Additionally, an excess emission is also evident in the profile extracted 
along the S-cavity and corresponds to the bright plume like feature shown by blue dotted circle 
in Figure |2d). To enhance the fluctuations seen in the regions of interest, we also plot the X-ray 
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photon distribution extracted from the rectangular strip region of the residual map (Figure |2}d) 
and is shown by red dashed curve in Figure |3] right panel). For comparison we also show the 
photon distribution extracted from similar region in the 2-d smooth model image (blue solid line). 
This figure clearly reveals a significant fraction of fluctuations in the X-ray photon count rates 
across the cavity positions in the residual image compared to the smooth distribution in the model 
image. An excess X-ray emission is also evident at about 40"and corresponds to the bright plum 
region in Figure Hid). 



4. Chandra Spectral Analysis 

A significant amount of variation, particularly in the central region, was evident in the surface 
brightness distribution of the Abell 1991 (Fig. H). With an objective to examine their relevance 
with temperature variation, we performed spectral analysis of the X-ray photons extracted from 
different regions of interest. For this, 0.5-7.0 keV X-ray spectra were extracted from each of the 
region using the CIAO task specextract and were grouped such as to get at least 20 counts 



in each bin. During t his study we fixed the 
= 2.46 x 10 20 cm~ 2 ) (IDickey and Lockman 



iy_drogen column density at the Galactic value (N H 



1990 ) and the redshift at z = 0.0587. The metal 



abundance and temperature were allowed to vary during the fit. 

To determine the global spectral properties of diffuse gas in Abell 1991, we extracted a 
combined spectrum from within 180"(~ 210kpc), excluding central 5"region and the resolved 
point sources. A total of 85,844 background-subtracted counts were extracted in the range between 
0.5— 7.0keV. Resultant spectrum was then fitted with a single temperature (phabs*apec) 
model which resulted in x 2 — 582 for 312 dof. Careful look at the residuals revealed that 
the single temperature model is not adequate to constrain the observed spectrum. Therefore, 
we tried with a double temperature model phabs (apec + apec) , which improved the fit 
to x 2 /dof = 372/241. The best fit lower and higher temperature values are found to be ~ 
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Fig. 1. — 0.5-7.0 keV raw Chandra ACIS-S3 image of Abell 1991, overlaid on which are the 21 
point sources detected within the chip S3. 



Table 1: Spectral properties of X-ray cavities, central knots, diffuse component (hot gas) and the 
central source 



Regions 


]\Tintr 


Best-fit 


kT 


kT 


Abund 


X 2 /dof 




(10 20 cm~ 2 ) 


Model 


(KeV) 


(KeV) 


z 




N-cavity 




(phabs*apec) 


i 7 7+0.19 
L -l ' -0.12 






47/62 


S-cavity 




(phabs*apec) 


i cro+0.05 
i - do -0.06 






178/137 


Excess emission reg. 




(phabs*apec) 


z - UD -0.07 




o.6i±S;g 


180/156 


Central source 


fi 9 c i +a02 

u - yj -0.03 


(phabs*phabs*apec+pow) 


n qq+o.03 

u - yo -0.02 




n 07+0.07 

u -°' -0.05 


70.20/75 


Total diff emission 




(phabs*apec+apec) 


±.JU_ 22 


6.84l 2 ;^ 


0.83+°i 


372/241 


Three knots 


7 sn +a02 


(phabs*phabs*apec) 


n 04+0.03 

u - y ^-0.03 




u - zd -0.04 


54/66 



Note Col. 1 - regions of interest used to extract the spectra in the energy band 0.5-7.0keV, Col. 2 - hydrogen column 
density intrinsic to each of the source region, Col. 3 - best-fit model, Col. 4 & 5 - low and high temperature values of 
the best fitted spectra, Col. 6 - metallicity values and Col. 7 - Goodness of fit. 
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Right ascension Right ascension 



Fig. 2. — (a) ACIS-S3 exposure corrected, background subtracted 0.5-3.0 keV Chandra image of 
Abell 1991. (b) Chandra 0.5-3.0 keV unsharp masked image of Abell 1991 derived after subtract- 
ing a 5o" wider Gaussian kernel smoothed image from that smoothed with a narrow 3a Gaussian 
kernel. This figure reveals substructures in the central region, (c) Elliptical 2-d beta model sub- 
tracted residual image exhibiting a pair of X-ray cavities and an excess X-ray emission region, (d) 
Same as that in (c) overlaid with VLA FIRST survey radio contours (white). The radio contours 
are at 0.5, 7.3, 14, 21, 27.23 mJy/beam levels, image rms is ~ 145 fiJy. Regions of interest are also 
shown in this figure. 
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Fig. 3. — (left panel): 0.5-5.0 keV azimuthally averaged surface brightness profile of the undis- 
turbed X-ray photon distribution (open triangles) within Abell 1991. Continuous line in it repre- 
sents the best fit /3-model. We neglected the regions of X-ray cavities as well as excess emission 
during this fit. Surface brightness profiles of the X-ray photons extracted along the conical sectors 
0-140° covering N-cavity (red hexagon) and 220-230° covering S-cavity (blue filled triangles) are 
also shown in this figure. The profile extracted along 220-230° conical sector also reveals an excess 
emission at ~ 40". (right panel): X-ray photon distribution within the 0.5-5.0 keV background 
subtracted 2d /3-model image (blue solid line) and those in the 2-d (3 model subtracted residual 
image (red dashed line). Significant depressions at the locations of both the cavities and an excess 
emission at ~ 20"-70"are clearly evident in this figure. 
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Fig. 4. — 0.5-7.0 keV radial temperature profile of the hot gas distribution in Abell 1991. This 
was obtained by fitting a single-temperature apec model to the spectra extracted form 1 1 different 
annuli centered on the X-ray peak of Abell 1991 cluster. Errorbars show the 90% confidence level. 

5 ' 

Knotl 




Fig. 5. — 0.5-2.0 keV exposure corrected background subtracted adaptively smoothed Chandra 
image of central 15"region of Abell 1991. This figure clearly revels three different knots in central 
region. Spectra representative from each of the knot have been extracted separately. 
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1. 961^22 keV and ~ 6. 84^3^4 keV, respectively, for the two c omponents and the c orresponding 



1+2.06 



Sharma et al 



(120041) fitting double 



X-ray luminosity is found to be 7.65 x 10 43 erg s" 1 (Table [U). 
temperature (MEKAL) model to the X-ray photons extracted from within 190"arrive at temperature 
values equal to 1.51 ±0.09 ke V and 7.87 ± 0.96 k eV. Our estimate of lower temperature value 



is higher than that deri ved by 



Sharma 



reasonable agreement. IStewart et al. 



et al 



(|2004f) . while the higher temperature value shows a 



19841) form deprojection analysis of total spectrum over 



Einstein 0.5-3.0 keV could find the temperature kT~1.6 keV. 



4.1. Azimuthally Averaged Spectral Analysis 



A projected radial temperature profile was also derived for Abell 1991 by analyzing spectra 
extracted from 1 1 different annuli centered on the cluster core. Widths of the annuli were 
adjusted such as to achieve roughly same S/N on each of the annulus. The background- subtracted 
0.5-7.0 keV spectrum extracted from each of the annuls was then fitted with an apec thermal 
plasma component modified by the intrinsic absorption (phabs*apec). We kept temperature 
and metallicity free during this fit. Profile of the best fitted temperature along with associated 
errors is shown in Fig. @] Like in other cooling flow galaxies, temperature profile derived for Abell 
1991 also shows a positive temperature gradient in the inner region, reaches to a maximum of 2.63 
keV at about ~7 6 kp c a nd then declines in the o utward direction, consistent with that reported by 



Sharma et al. 



(2004) and 



Cavagnolo et al 



(2009) 



4.2. Cavities and excess emission region 

To perform a more detailed analysis of the temperature structure of the hot gas corresponding 
to the cavities and the excess X-ray emission regions, we extracted spectra representative of both 
the cavities (shown by yellow ellipses in Fig. [2]d) as well as from the excess X-ray emission 
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region (shown by dotted blue region in Fig. |2}1). 0.5-7.0 keV spectra from these regions were 
then fitted independently with the single temperature model (phabs*apec) within XSPEC, with 
the column density set to the Galactic value and keeping temperature and abundance values free. 
This yielded the best fit temperature values for N and S cavities to be equal to l-77;to;]J keV and 
l-53lo!o6 keV, respectively, while that for the excess emission region was found to be 2.06l ; 7 
keV. The metal abundance associated with the excess emission is found to be 0.6l7_o!o9 Z©. This 
means the region of excess emission is hotter compared to the cavity regions, which themselves 
are relatively hotter than the undisturbed gas surrounding the cavities. It is likely that the X-ray 
photons extracted from cavities may also include photons from the foreground and/or background 
of the cavities. However, our projected and deprojected analysis of the X-ray photons from this 
system, particularly in the cavity locations, does not show significant variations in the estimated 
values of temperature as well as luminosity. Therefore, we believe that the effect of such emission, 
if any, would be negligible. 



4.3. Central Source 

To examine spectral properties of X-ray emission from the nuclear region of Abell 1991, 
we extracted a 0.5-7.0 keV spectrum from the central 5 "region centered on the X-ray peak of the 
cluster. This spectrum was fitted with a single temperature apec model, however, it resulted in 
poor fit particularly at higher energy range. Therefore, we added a power-law component to this 
model. The best-fit temperature obtained from the double component model is 0. 931^02 keV and 
is much cooler than the surrounding ICM, while the abundance is found to be equal to 0.37^57 
Z . The power-law index in the present case is T=1.4, indicative of the sufficient hard structure 
of the central source. The unabsorbed 0.5-7.0 keV flux of the central source was found to be 
5.13 x 10~ 13 erg s" 1 cm" 2 leading to a luminosity of 3.95 x 10 42 erg s -1 . The best fitted spectral 
parameters of the central source are given in Table [T] 
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4.4. X-ray knots 



Exposure corrected, background subtracted adaptively smoothed 0.5-2.0 keV Chandra image 
revealed three different knot-like features in the central 15 "region of the Abell 1991 tha t are ab out 



Sharma et al. 



(2004). 



10"off the X-ray peak (Fig. [5]) and are consistent with those reported by 
To examine properties of these knots, we extracted a combined spectrum from these knots and 
exported it to XSPEC for spectral fitting. We fitted the spectrum with an apec model including 
an additional absorption component. The best fit resulted in the tem perature equal to 0.9 4±0.03 



Sharma et al. 



(120041) with 



keV and metallicity equal of 0.251q:o4 Z (Table[B, while those by 
(WABS+MEKAL) are 0.83±0.02 keV and 0.46i°;o| Z , respectively. Our estimates are slightly 



higher than those of 



Sharma et al. 



(|2004l) . Spectra representative from each of the three knots 
were also extracted and fitted separately. Spectra of knots 1 and 2 were constrained with simple 
(apec) component, however, that for knot 3 required an addition absorption component. 



4.5. X-ray Point Sources 



The wavelet based detection algorithm wavelet ect within (CIAO) enabled us to detect a 
total of 21 discrete sources within ACIS-S3 chip of 0.5 -7.0 keV Chandra image of Abell 1991. 
As flux of the X-ray photons from individual source was not enough to perform spectroscopy 
of individual source, therefore, we de rived X-ray colors of these sources with an objective to 



investigate their structural properties (|Prestwich et al 



2009; 



Vagshette et al. 



2012). For this we 



extracted X-ray counts for each of the source in three different energy bands: soft (S, 0.5-1.0keV); 
medium (M, 1.0-2.0keV) and hard (H, 2.0-7.0 keV) using the task dmextract. Then the 
hardness ratios of the individual source were derived using definitions H21=(M-S)/(M+S) and 
H31=(H-S)/(H+S) and are listed in Tabled We then plotted the X-ray color-color plot for the 
resolved sources and is shown in Figure |6] where are the sources are found to lie along the 
diagonal region covering hardness ratios in the range between [H21,H31] = [-0.55, -0.52] to 
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Table 2: Properties of the point sources detected within ACIS-S3 chip 


No. 


RA 


DEC 


Area 


H21 


H31 Count Rate 


Luminosity 








(arcsecr) 




(io- 


3 counts s x ) 


(10 40 erg s" 1 ) 


1 


14:54:41.017 


+18:37:29.06 


34.86 


0.47 


0.20 


0.548 


3.51 


2 


14:54:39.282 


+18:35:24.72 


24.21 


0.65 


0.58 


0.992 


6.35 


3 


14:54:39.137 


+18:37:54.61 


30.98 


0.27 


0.60 


0.705 


4.51 


4 


14:54:37.835 


+18:39:55.36 


28.08 


0.67 


0.71 


0.627 


4.01 


5 


14:54:35.101 


+18:41:03.35 


15.49 


-0.12 


-0.55 


2.611 


1.67 


6 


14:54:32.035 


+18:39:59.72 


14.52 


0.56 


0.59 


1.305 


8.35 


7 


14:54:31.704 


+18:39:47.34 


13.56 


0.04 


-0.42 


1.436 


9.19 


8 


14:54:25.296 


+18:36:18.49 


20.33 


0.31 


-0.25 


0.913 


5.85 


9 


14:54:22.154 


+18:41:04.92 


19.37 


0.16 


-0.03 


1.384 


8.86 


10 


14:54:42.809 


+18:39:02.64 


38.73 


-0.29 


-0.29 


0.600 


3.84 


11 


14:54:36.171 


+18:41:18.91 


99.73 


-0.52 


-0.55 


1.593 


1.02 


12 


14:54:34.645 


+18:37:26.32 


25.17 


0.07 


-0.04 


1.1230 


7.19 


13 


14:54:32.214 


+18:40:52.96 


56.16 


-0.10 


-0.13 


1.149 


7.35 


14 


14:54:23.975 


+18:35:23.19 


77.46 


0.00 


-0.04 


0.913 


5.85 


15 


14:54:23.306 


+ 18:35:04.75 


25.17 


0.12 


-0.29 


0.809 


5.18 


16 


14:54: 1 /.959 


. lo.'l c . A^l AA 

+ 18:35:02.00 


C\A OA 

94.89 


A A A 

0.49 


0.65 


2.454 


1 Cl 

1.57 


17 


14:54:16.001 


+18:39:57.66 


44.54 


0.06 


-0.33 


1.514 


9.69 


18 


14:54:29.249 


+18:34:42.54 


80.37 


0.27 


0.38 


1.044 


6.68 


19 


14:54:23.751 


+18:38:10.14 


39.70 


0.10 


0.10 


0.809 


5.18 


20 


14:54:07.675 


+18:36:42.88 


29.05 


0.18 


-0.08 


0.600 


3.84 


21 


14:54:06.714 


+18:36:59.70 


84.24 


-0.25 


0.29 


0.887 


5.68 



Notes: Col. 2 & 3 - position coordinates of the resolved sources, Col. 4 - area information for each of the source, the 
hardness ratios H21 and H31 are given in the Col. 5 & 6, respectively. Col. 7 & 8- count rates and X-ray luminosities 
(Lx ) of individual source, respectively. 
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[+0.68, +0.71]. Among them, four are harde r sources with hardn ess ratios > [+0.50, +0.50] and 



probably represent strongly absorbed AGNs (jRandall et al. 



2004) 



5. Discussion 



5.1. Cavity Energetics 



Comparison of the cavity power with the X-ray luminosity of the ICM within cooling radius 



allow us to examine balance between the radiative 



(Birzan et al. 


2004; 


Dunn et al. 


2005; 



Raffertv et al. 



oss fr om the ICM and the cavity heating 



200 6). Therefore, an attempt was made to 



estimate power of both the cavities b y measuring pdV work done b y the jets in inflating these 
cavities and by estimating their ages (|McNamara and Nulsenl 120071) . To investigate properties 
of X-ray cavities and to quantify energetics of the AGN outburst it is required to know the ICM 
density, temperature and pressure at the cavity location. For this we used azimuthally averaged 
temperature profile of the X-ray photons derived above in §4. 1 . Energetics of both the cavities 



were then determined following standar d method and using the relation E, 



(Birzan et al 



2004; 



Raffertv et al 



cavity 



71 



71- 



20061) . where p is thermal pressure of the ICM at the cavity 
positions and V is the volume of each cavity. For the relativistic fluid using 71 = 4/3 this 
equation becomes i^vity = 4pV. Volumes of both the cavities were estimated assuming that 
they are symmetric about the plane of the sky with centers lying in the plane perpendicular to the 
line-of-sight and passing through the central AGN. Then the expression for volume of the cavity 
becomes V = AnabR\ J3, where a and b are the projected semi-major and semi-minor axes of 
the cavities, respectively, and Ri os is the line-of-sight distance between the nuclear source and the 
center of each cavity. The measured dimensions of both the cavities along with their volumes are 
given in Table |3j 



We adopted all the three methods proposed by 



Birzan et al. 



(|2004r) to estimate ages of the 
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H21=M-S/M+S 

Fig. 6. — X-ray color-color plot of all the 21 point sources detected within 0.5—7.0 keV ACIS-S3 
chip. 



Table 3: Properties of the X-ray cavities within Abell 1991 



Cavity Parameters 


Ncavity 


Scavity 


a (kpc) 


16.8 


13.3 


b (kpc) 


5.5 


6.2 


R (kpc) 


12.4 


11.5 


t Cs (yr) 


1.8xl0 7 


1.8xl0 7 




3.8xl0 7 


2.9xl0 7 


trefilliy*) 


6.8 xlO 7 


5.9xl0 7 


Pcavity (d"g S ) 


3.8 xlO 43 


4.8 xlO 43 



Note : a , b - semi major and semi minor axis lengths of the cavities, R projected distance from the center of the 
cluster. 
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cavities, namely, (1) the cavity reaches the projected distance R from the center of the cluster 
in sound crossing time t s = R/c s = Rj y/^kT/ ^,m p ; where 7=4/3, k is Boltzmann constant, 
T is temperature, fi is mean molecular weight and m p is mass of the proton; (2) cavities are 
rising buoyantly to the present position in time t buoy = Rj ^j2gV/SC, where g = GM/R 2 is the 
gravitational acceleration at the cavity position R, V is volume of the cavity, S is the cross-section 
of the cavity, and C = 0.75 is the drag coefficient; and (3) time required for the gas to refill the 



displaced volume of the cavity as it rises in the time given by t re f ~ 2y/r/g, where r is radius of 
the cavity. 

The buoyancy rise time is the time taken by the cavity to reach its terminal velocity that 
depend on the drag forces and in the case of clearly detached cavities provides better estimate of 
the age. The sound crossing time is the time taken by the bubble to travel the projected distance 
from central AGN to its present position at the sound speed, while the refill time is the time taken 



by a bubble to rise buoyantly through its own diameter starting fr om the rest 



Vpically, the ages 



2004), with the 



estimated by all the three assumptions agree within a factor of 2 (B irzan et al. 
buoyancy time lying between the sound crossing time and the refill time. We estimated ages of 
both the cavities using all the three methods discussed above and are given in Table |3] Though 
all the three estimates do not vary significantly, which among these represent the best estimate of 
age is still not clear. Therefore, in the present case we took average of all the three estimates and 
derived the cavity power using relation P cav i ty = <^ aYlt > . Resultant values of total powers for the N 
and S-cavities in the present study are found to be equal to 3.8 x 10 43 erg s~ x and 4.8 x 10 43 erg s -1 , 
respectively, and are given in Table |3] As these estimates do not include possible hydrodynamical 
shocks, therefore, provides the lower limit to the total power of the AGN. 
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Fig. 7. — (left pane l): Cavity heat i ng po wer vs the bolometric X-ray luminosity within cooling 



radius adopted from 



Raffertv et al. 



(|2006l) . The diagonal lines denote cases where heating power 
Licm) fox pV, ApV and \6pV total enthalpy levels of the cavities 



equal cooling power (P cav 
from top to bottom. For comparison, we show a balance between the two for the case of Abell 
1991 (open red triangle), (right panel): Comparison of th e balance between cavit y heating power 



Cavagnolo et al 



d2010h . The red long 



(P C a«) v s 1.4 GHz radio power for the sample studied by 
dashed line represents the best fit relation for their sample of giant ellipticals (gEs), while the short 
dashed blue lines represent the upper and low er limits for their be st fit relations. For comparison, 



we also show the best fit relation obtained by 



Birzan et al. 



(120041) for galaxy c l usters (solid black 



line). Abell 1991 occupies position little above the best-fit line by 



Birzan et al 



(2004) 
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5.2. Cavity Scaling Relations 



The re are mounting evidences stating that the central dominant gal axies (CDGs) host radio 



sources (Birzan et al. 



2004; 



Raffertv et al. 



2006 



Cavagnolo et al. 



2010|) and exhibit profound 



effect on the X-ray emitting ICM surrounding such systems. The ambipolar jets originating 
from the central source injects energy in to the ICM at regular intervals and hence results in the 
formation of under dense bubbles on the opposite sides of the nucleus which rises buoyantly 
in the cluster and expand until they reach pressure equilibrium with the surrounding ICM 



(Blandford and Rees 



19741) . But whether the energy deposited by the central AGN is sufficient to 



account for the buoyant rise of the bubbles in the ICM is not clear. 

With an objective to investigate association of this cooling flow galaxy with the central 
radio source, we have made use of the radio data available in the archive of the Very Large 
Array Faint Images of the Radio Sky at Twenty-Centimeters (VLA FIRST) survey. The radio 
emission contours derived from the analysis of radio data are overlaid on the unsharp masked 
X-ray image of Abell 1991 and is shown in Figure |2](d). This figure confirms the coincidence of a 
po int-like source with t he central dominant galaxy NGC 5778 and is consistent with that reported 



by 



Sharma et al. 



(120041) . however, shows an offset of 10"relative to the X-ray peak of the cluster. 
Though this radio source shows an association with the S-cavity, no such association is evident 
for the case of N-cavity. Such cavities with n o apparent spatial coincidence with the 1 .4 GHz 



2000; 



McNamara et al. 



2001 ) and are 



radio contours are referred as "ghost cavities" (jFabian et al 
thought to be formed by the previous outbursts. 

To investigate whether the energy deposited into the ICM from the central radio source is 
sufficient to account for the heating, we compare the mechanical power of the AGN with the 
cooling power of the ICM. AGN power, as discussed above, was estimated as the pdV work 
required to inflate a cavity against the surrounding pressure of the hot ICM. For the case of 
Abell 1991, the total output power of the radio source derived from the measured values of 
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pressure and volume of the X-ray cavities is found to be ~ 8.64 x 10 43 erg s~\ while the X-ray 
luminosity within the cooling radius (~ 117 kpc) is found to be L IC m(< r cooz)~ 6.04xl0 43 erg 
s _1 . Therefore, in this case, on the basis of simple energy arguments, it is evident that the radio 
source hosted by Abell 1991 on average is capable to deposit enough energy into the ICM so 
as to offset the cooling g as. We further check this balance by comparing its relevance with that 



of the sample studied by 



Raffertv et al. 



(|2006l) . For this we make use of their plot between the 
cavity power of the central AGN against total radiative luminosity of the intracluster gas within 
the cooling radius (Figure [7] left panel). The thick diagonal lines in this figure represents the 
P cav = Ljcm for the total enthalpy levels equal to pV, 4pV and 16pV of the cavities. In this plot 
we indicate the Abell 1991 by red open triangle and is found to occupy a position on the ApV line. 
This in turn imply that the total power available within the cavities of Abell 1991 is sufficient to 
balance the entire radiative loss of the ICM within the cooling radius. 



Cavagnolo et al 



(120101) observed a strong power-law relationship between the mechanical 
power stored within the cavities and the 1 .4 GHz radio power of the central source. We made 
use of this relationship to compare 1 .4 GHz radio power {P\ a ght) and cavity power (P m „) of 



the Abell 1991 using the samples studied by 



Cavagnolo et al. 



( 20101) and 



Birzan et al. 



( 20041) and 



is shown in Figure [7] (r ight panel). The red long dashed line represent the best fit relation of 



Cavagnolo et al. 



(120101) for the sample of giant ellipticals (gEs), while small dashed blue lines 



de notes the upper and lower limits for their best fit relations. We also show the best-fit relation 



of 



Birzan et al 



(|2004|) for a sample of galaxy clusters studied by them and is shown by thick 
continuous line. In this plot Abell 1991 occupies a position confirming the correlation between 
the two and hence imply that the radio source hosted by this system has a capability to reverse the 
cooling flow. We have also compared the jet power (Pj et ) of radio source host ed by Abell 1991 



with th at of the observed radio power (V radio) using the r elation presented by 



(|2010j) and found that, like other FR-I syste ms studied b y 



exhibit similar power capabilities. Further, 



Birzan et al 



Croston et al 



Cavagnolo et al. 



d2008|), Abell 1991 also 



(|2012l) have demonstrated that every 
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massive cluster hosts a radio source with its radio luminosity greater than 2.5 x 10 erg s Hz 
Compared to this, the radio source associated with Abell 1991 have 1.4 GHz luminosity equal to 
~ 3.02 x 10 30 erg s _1 Hz -1 , and hence support the idea of self-regulating radio-mode outburst. 
Thus, all the attempts made here to relate the jet/cavity power with the X-ray /radio luminosity 
imply that the AGN outburts from the Abell 1991 are comparable to those needed for quenching 
the cooling flow. However, deep, high-resolution, multi-frequency radio observations of this 
system are called for arriving at proper conclusion. 

6. Summary 

In this paper we have presented results based on the re-analysis of the Chandra X-ray 
observations of the galaxy cluster Abell 1991 with an objective to investigate properties of X-ray 
cavities and the interplay with the central radio source. Summary of the main results is: 

1 . Unsharp masked image as well as 2-d (3 model subtracted residual image reveals a pair of 
ellipsoidal shaped X-ray cavities at a projected distance of about 12kpc from the central 
region of Abell 1991 and are consistent with those seen in several other cooling flow 
galaxies like NGC 6338. In addition to these prominent cavities, an excess X-ray emission 
region is also evident on the southern side of the S-cavity. Surface brightness profiles 
derived for this system revealed depressions in the X-ray emission at a distance of about 
12kpc, confirming the presence of X-ray cavities. 

2. Spectral analysis of X-ray photons extracted from both the cavities gives temperature of 
gas corresponding to the N and S cavities equal to 1.77 ± 0.15 keV and 1.53±0.06 keV, 
respectively, while that of the excess X-ray emission region is 2.06± 0.09 keV. 

3. Like in other cooling flow clusters, temperature profile derived for Abell 1991 exhibits a 
positive temperature gradient, reaching to a maximum of 2.63 keV at ~ 76 kpc and then 
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declines in outward direction. 

4. 0.5—2.0 keV X-ray image of Abell 1991 reveals three different knot like structures in the 
central 15"region that are about 10"off the X-ray peak of the cluster. However, an additional 
absorption component was required to constrain combined spectrum of the knots. 

5. Total mechanical power of the cavity system in Abell 1991 is P ca ^ y =8.64x 10 43 erg s -1 , 
while the luminosity of the gas within cooling radius (r coo /=l 17 kpc) is L coo / ~ 6.04 x 10 43 
erg s" 1 . This is turn imply that the mechanical power injected by the central AGN through 
outbursts is sufficient to balance the radiative loss. 

6. Radio emission map derived using VLA FIRST survey data shows an apparent coincidence 
with the S-cavity and exhibit an offset of 10"relative to the X-ray peak of the Abell 1991. 

7. AGN heating, as traced by the power content of the X-ray cavities, is capable of balancing 
the radiative losses of the ICM in this system. However, high resolution, low-frequency 
GMRT observations of Abell 1991 are required to examine the association of X-ray cavities 
with the radio-jets. 

8. We have detected a total of 21 point-like sources within the chip 7, majority of which are 
like 1.4 Af accreting neutron star binaries while few appears as heavily absorbed harder 
sources. 
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